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ABSTRACT 

Narrow planetary rings are eccentric and inclined. Particles within a given 
ring must therefore share the same pericenter and node. We solve for the three- 
dimensional geometries and mass distributions that enable the Uranian a and 
(3 rings, and the Saturnian Maxwell and Colombo (Titan) rings, to maintain 
simultaneous apsidal and nodal lock. Ring self-gravity, interparticle collisions, 
and the quadrupole field of the host planet balance each other to achieve this 
equilibrium. We prove that such an equilibrium is linearly stable. Predictions for 
the Saturnian ringlets to be tested by the Cassini spacecraft include: (1) ringlet 
masses are of order a few x 10 19 g, (2) surface mass densities should increase 
from ring midline to ring edges, and (3) rings are vertically warped such that the 
fractional variation of inclination across the ring is of order 10%. Analogous pre- 
dictions are made for the Uranian rings. Simultaneous apsidal and nodal locking 
forces the narrowest portion of the ring — its "pinch," where self-gravitational and 
collisional forces are strongest — to circulate relative to the node, and introduces 
previously unrecognized time-varying forces perpendicular to the planet's equa- 
tor plane. We speculate that such periodic stressing might drive kilometer-scale 
bending waves at a frequency twice that of apsidal precession; such flexing might 
be observed over a few weeks by Cassini. 

Subject headings: celestial mechanics — planets and satellites: individual (Uranus, 
a and (3 rings, Saturn, Maxwell and Colombo rings) 

1. INTRODUCTION 

The a and (3 rings of Uranus are eccentric and inclined (Elliot et al. 1984; French et 
1991). The mean eccentricity, e, of each ring is (0.761 ± 0.004) x 10~ 3 and (0.442 ± 
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0.003) x 10~ 3 , respectively. The mean inclination, J, of each ring with respect to the equator 
plane of Uranus is 0.0152 ± 0.0006 deg and 0.0051 ±0.0006 deg, respectively. To maintain 
an observed eccentricity and inclination, each ring must be composed of particles that share 
the same longitude of pericenter and same longitude of ascending node. The alignments of 
apsides and of nodes must simultaneously follow from a delicate balance of forces due to the 
gravitational field of the planet, ring self-gravity, and interparticle collisions. Obtaining a 
complete equilibrium solution for a ring's three-dimensional geometry and mass distribution 
has been a goal of dynamicists for decades. The solution bears directly on ring ages and 
origins (see, e.g., Goldreich & Tremaine 1979a; Chiang 2003). 

The Maxwell and Colombo ringlets are Saturnian analogs of the Uranian a and (3 rings 
(Porco 1990). 1 The mean eccentricities of the Saturnian ringlets, but not their mean incli- 
nations, are measured from Voyager spacecraft observations. Table 1 summarizes observed 
ring parameters. 

Shepherd satellites are expected to confine narrow rings and maintain the latter's sharp 
edges (Borderies, Goldreich, & Tremaine 1982), but no satellite has yet been observed to 
accompany any of the aforementioned ringlets. This state of affairs promises to change with 
the arrival of the Cassini spacecraft to Saturn in 2004. Narrow rings and their attendant 
shepherds furnish the most accessible laboratories we have for studying disk-satellite interac- 
tions. Their study informs problems on grander scales, including migration of planets within 
circumstellar disks (see, e.g., Chiang 2003). 

The steady-state condition of apsidal alignment, combined with measurements of the 
ring eccentricity profile, e(a), where a is the semi-major axis of a ring streamline, can be 
employed to calculate the ring surface density, S(a). First attempts at this calculation 
account for the planetary quadrupole field and ring self-gravity, but omit the effects of 
interparticle collisions (Goldreich & Tremaine 1979b, hereafter, the "standard self-gravity 
model"). Chiang & Goldreich (2000, hereafter CG00) restore the latter effects to demonstrate 
how collisions near ring boundaries can dramatically raise the ring masses required for apsidal 
alignment. Surface densities derived by accounting for pressure forces exceed those derived 
without pressure forces by 1-2 orders of magnitude; the larger masses can be reconciled 
with occultation measurements and with classical theories of ring shepherding (CG00). A 
more careful treatment of pressure forces in the ring interior is undertaken by Mosqueira & 
Estrada (2002), who reach the same qualitative conclusions as those of CG00. 

All the works cited above consider only apsidal alignment and neglect nodal alignment. 
Nodal alignment, but not apsidal alignment, is studied by Borderies, Goldreich, & Tremaine 



The Colombo ringlet is also known informally as the Titan ringlet. 
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(1983a), who consider ring self-gravity and planetary gravity, but neglect collisions. Ring 
eccentricities are set to zero in their analysis. 

The ring's true surface density profile must be simultaneously reconciliable with both 
the alignment of apsides and the alignment of nodes; the horizontal structure of a narrow 
ring is entwined with its vertical structure. This paper seeks to simultaneously treat apsidal 
and nodal alignment while accounting for the full panoply of forces due to the planetary 
quadrupole field, ring self-gravity, and interparticle collisions. In §2, we derive equilibrium 
ring surface densities and vertical geometries that lock the apsides and nodes of a given 
ring. We apply our solutions to the a and (3 rings of Uranus, and the Maxwell and Colombo 
ringlets of Saturn. In §3, we present a proof that circular, nodally locked rings are linearly 
stable to perturbations to their inclinations and nodes. The beginnings of such a proof can 
be found in Borderies, Goldreich, & Tremaine (1983b); here, we state the arguments more 
completely and explicitly. In §4, we discuss our results, highlighting the future impact of the 
Cassini spacecraft on studies of narrow rings and unresolved theoretical issues. 



2. EQUILIBRIUM 

Our procedure for deriving the mass and 3-dimensional structure of a narrow ring is 
summarized as follows. The range of semi-major axes spanned by the ring, the eccentricity 
profile [e(a)], and the mean inclination (J) are assumed to be given. From e(a), we com- 
pute the surface density profile, S(a), by enforcing apsidal alignment across the ring and by 
accounting for planetary oblateness, ring self-gravity, and interparticle collisions. 2 This com- 
putation is described in detail in §2.1. Next, from S(a) and /, we compute the inclination 
profile, 1(a), by enforcing nodal alignment across the ring and by accounting for planetary 
oblateness and ring self-gravity but not interparticle collisions. The computation of 1(a) is 
described in §2.2. Finally, in §2.3, we gauge a posteriori the extent to which our neglect of 
collisional stresses in the derivation of 1(a) is justified. 

Our computational procedure assumes the ring is characterized by a time- independent 
argument of pericenter, u. This means the narrowest portion of the ring (its pericenter) 
lies at a fixed angle relative to the ring's node on the equator plane of the planet. Strictly 



2 The Colombo ringlet has the added complication that it inhabits a 1:0 apsidal resonance established by 
the Saturnian satellite, Titan; the ring's apsidal precession frequency matches the mean motion of Titan. It 
is straightforward to show that the contribution of Titan to the differential apsidal precession rate across 
the Colombo ringlet is small compared to the differential apsidal precession induced by the quadrupole field 
of Saturn. We will therefore neglect perturbations by Titan on the Colombo ringlet in this paper. 
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speaking, the assumption of stationary uo is invalid, since 

Co = Q - (l = 2u ^ 0. (1) 

The common precession rate of the longitude of pericenter of all apsidally aligned streamlines, 
oj, is set by the quadrupole field of the central planet and is positive (in the direction of 
increasing true anomaly). The common precession rate of the longitude of ascending node of 
all nodally aligned streamlines, Cl, is also set by the planetary quadrupole field, is identical 
in magnitude to oj, and is negative. That Co ^ implies that the pericenter of the ring 
circulates relative to the node on the planet's equator plane; since self-gravitational forces 
are largest at pericenter, the contribution to differential nodal precession from ring self- 
gravity cannot be time-independent. Our hope is that this complication does not introduce 
serious errors into our calculation of the time-averaged ring geometry; we point out below 
where our assumption of fixed argument of pericenter is employed, explore in a preliminary 
but quantitative manner the consequences of its violation in the Appendix, and speculate 
about its true effects in the discussion section (§4). 



2.1. Surface Density 

We divide the ring into a set of 2N wires and compute the mass of each wire by imposing 
the steady-state condition of apsidal alignment. The 2N wires are equally spaced in semi- 
major axis; N wires lie at semi-major axes smaller than that of the ring midline, and N 
wires lie at greater semi-major axes. Our calculation differs from that of CG00 only in the 
treatment of collisional forces. We restrict ourselves here to describing this difference; for 
background details, see CG00. 

CG00 highlight the importance of interparticle collisions at ring boundaries and account 
crudely for their effects by introducing force terms in the equations of motion for "endwires" 
located within a collisional mean free path, A, of either ring boundary. No account is made 
of collisional stresses in the ring interior (regions displaced many mean free paths from either 
boundary) in the quantitative model they present. Here we improve upon their calculation 
by making such an account, although still in a crude and prescriptive manner. We take the 
solution of CG00 for S(a) and derive from it a collisional acceleration everywhere within the 
ring: 

C( , /H -!^_^_ f , (2) 
L da 1 — q e cos / 
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where f points in the radial direction, q e = ade/da is the dimensionless eccentricity gradient, 
and / measures true anomaly. The factor 1 — q e cos / is proportional to the local radial 
spacing between wires (see, e.g., Goldreich & Tremaine 1979b). Following CGOO, we adopt 
a constant q e (a) = q e across the entire ring. The velocity dispersion, c, and a dimensionless 
factor, S, are described as follows. 

We prescribe the following profile for the square of the velocity dispersion, 



where q, q,, and w r > are the same constants as employed by CGOO, and aj, is the 
semi-major axis of the ring boundary that is nearest a given wire at a. According to (3), 
the velocity dispersion decreases from either ring boundary towards the ring interior over 
a lengthscale w r . The physics underlying the enhancement of velocity dispersion near ring 
boundaries is elucidated by Borderies, Goldreich, & Tremaine (1982) and estimates for q, 
c fe , and w r are derived by CGOO. The actual values we employ are contained in Table 2. 

The function 



is a softening parameter that we impose because the true acceleration due to particle collisions 
near a ring boundary is likely overestimated by the usual hydrodynamic expression for the 
acceleration due to pressure gradients [— V(Sc 2 )]. Within a few mean free paths, ~A, 
of the ring edge, the hydrodynamic approximation breaks down and particles behave more 
ballistically with less regard for large-scale gradients in the surface density. The parameter 
S quickly grows from e -2 to unity as we recede from the ring edge towards the ring midline. 

The collisional acceleration, C, is inserted into Gauss's equation for a given ring par- 
ticle's apsidal precession rate, du/dt, and averaged over true anomaly to yield, for a given 
wire, {du/dt)c- This collisional contribution to the wire's precession rate adds to other con- 
tributions due to planetary oblateness and inter-wire gravity. Expressions for the latter two 
contributions are supplied by CGOO. The condition of apsidal alignment requires that pre- 
cession rates of all wires be equal to the precession rate of a test particle at the ring midline, 
and yields 2N linear equations for the 2N wire masses. The solution of this linear system 
generates a new surface density profile, S'(a). The entire calculation is then repeated: from 
this new surface density, we calculate a new collisional acceleration profile [C'(a, /)], a new 
set of wire precession rates due to pressure forces [{du>/dt)c], and a new set of 2N linear 
equations. In this way, the solution for the surface density is iterated until convergence is 




(3) 



S = exp[-2A/(|a-a 6 | + A)] 



(4) 
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achieved. 

In practice, for model parameters appropriate to the Uranian a and f3 rings and the 
Saturnian Maxwell and Colombo ringlets (see Table 2), we find that the solution of 2N = 
2000-6000 wire masses converges after ~10 iterations. Because the solution is reflection- 
symmetric about the ring midline, only N = 1000-3000 linearly independent equations need 
be solved for N distinct wire masses. The derivative in equation (2) is computed numerically 
using a Savitzky-Golay smoothing filter of order 2 and having a width of 20 wires (Press 
et al. 1992). The system of 2N linear equations is solved using subroutine DGESV of the 
LAPACK (Linear Algebra Package) software library. Because many of our parameters such 
as c and S are only order-of-magnitude estimates, our profiles are probably accurate to 
factors of a few, at best. Nevertheless, by accounting not only for interparticle collisions 
near ring boundaries but also for collisions within the ring interior, we may explore the 
qualitative effects of incorporating the latter and make a first-cut correction to the solutions 
obtained by CG00. 

Computed surface density profiles for the a, (3, Maxwell, and Colombo ringlets are 
displayed in Figures 1-4, respectively. Near the ring midline, our solutions require less mass 
than those of CG00, a consequence of local pressure gradients that compress the ring and abet 
ring self-gravity. By contrast, our computed peaks in surface density near ring boundaries 
are larger than those of CG00. The reason for this is as follows. As a given peak in surface 
density is approached from the ring midline, steep and inward-directed pressure forces must 
be balanced by the outward gravitational attraction of massive endwires on the far side of 
the peak. These inward-directed forces are neglected by CG00; accounting for them here 
leads to endwire masses larger than those obtained by CG00. 

These same qualitative conclusions are reached by Mosqueira & Estrada (2002), who 
employ various prescriptions for collisional stresses that differ from ours. Their surface 
density profiles deviate from ours by large amounts, sometimes by more than an order 
of magnitude, reflecting the sensitivity of the shape of the profile to choice of boundary 
conditions. Some consolation may be had in the finding of Mosqueira & Estrada (2002) 
that the total ring mass, M, is less sensitive to this choice. The chief shortcoming of all 
calculations of the surface density profile, including our own, is the prescriptive and slightly 
arbitrary description of the ring boundary. We have assumed here that q e is constant across 
the entire ring, even in boundary regions resonantly perturbed by shepherd satellites where 
streamlines can no longer be described as simple ellipses. See CG00 for a discussion of this 
point. 
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Fig. 1. — Computed surface density profiles for the a ring of Uranus, for q e (a) = q e = 0.472. 
Including pressure forces in the ring interior lowers the surface density near the ring midline 
and raises it near the ring boundary, as compared to the model by CG00 who neglect interior 
pressure gradients. The total mass of the ring, M, hardly changes between models, however. 
The bottom dashed line corresponds to the standard self-gravity (SSG) model of Goldreich & 
Tremaine (1979b); it predicts surface densities too low to accord with observations (CG00). 
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Fig. 2. — Same as Figure 1, but for the j3 ring of Uranus (q e = q e = 0.370). 
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Fig. 3. — Same as Figure 1, but for the Maxwell ringlet of Saturn (q e = q e = 0.46). 
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Fig. 4. — Same as Figure 1, but for the Colombo (Titan) ringlet of Saturn (q e = q e = 0.44). 



2.2. Inclination Profile 



Given I and S(a), we compute /(a) by imposing the steady-state condition of nodal 
alignment. We now discretize the ring into 2N — 1 wires that share the same longitude of 
ascending node and solve for the orbital inclination of each wire with respect to the equator 
plane of the planet. Note that whereas the calculation of the surface density in §2.1 employs 
an even number of wires (2iV), here we employ an odd number for the calculation of the 
inclination profile. The reason for this change of procedure is to have the iV th wire lie on the 
ring midline and to assign In = I- With this reference point defined, the inclinations of the 
N — 1 wires that lie to one side of the midwire may then be computed. Inclinations of the 
N — 1 wires that lie on the other side of the midwire follow from reflection anti-symmetry 
about the midwire. 

The nodal precession rate of the j th wire relative to that of a particle on the ring midline 
and due to the quadrupole field of the central planet is given by 



where a and n are the semi-major axis and mean motion, respectively, of a particle on the 
ring midline, J 2 is the usual dimensionless measure of the planetary quadrupole field, R p is 
the planetary radius, and () denotes a time-average over one orbit. Here j runs from 1 (the 
wire closest to the planet) to 2N — 1 (the wire furthest from the planet). 

The corresponding differential nodal precession rate due to inter-wire gravity is 



2JV-1 

TrMpI ^ a k - a,j 



where 



_ h — Ij ,„s 
qkj = qjk = a , (7) 

0,k — dj 

M P is the mass of the planet, m k is the mass of the k th wire, and B is a dimensionless 
function of order unity that depends on the eccentricity profile of the ring, the inclination 
profile (qjk), and the arguments of pericenters of the various wires. Following our assumption 
of a time-independent oo (see the discussion preceding §2.1), we will take B = 0.77 constant. 



Though B actually varies with time because Co ^ 0, its maximum range of variation is less 
than a factor of 2, as we derive in the Appendix. 

The N distinct wire masses computed in §2.1 furnish the set of {rrij} used in equation 
(6) from j — 1 to iV; the remaining N — 1 wire masses are assigned by reflection symmetry 
across the ring midline. 

Steady nodal alignment implies 



A Q (Q) 3 + A G (fi) j = {j = 1, . . . , N - 1} , (8) 

where we have neglected the contribution from interparticle collisions; the validity of this 
omission is tested in §2.3. Embedded in equation (8) are N — 1 independent equations for 
the N — 1 distinct values of {Ij}. 

We proceed to massage equation (8) into a form that permits easy solution. A dimen- 
sionless form of the j th sub-equation reads 



2JV-1 

— hi 



J-N AH(2N-1)M fgy f ay rLk 



2N-1 2l7rM P iJ 2 \R P J \AaJ k - j 

where Aa = a27v-i — a\ is the total width of the ring near quadrature, and hk = mk/M is the 
fraction of the ring's total mass contained in wire k. Now each q^j is a linear combination 
of the N — 1 independent variables, 

Xi = q iii+ i = a— {i = 1, . . ., N - 1}, (10) 

as in 



j+i-i 

Qjj+i-m E x i- ( n ) 

i=j 

Furthermore, for ring self-gravity to balance the planetary quadrupole field, the inclination 
profile must be reflection anti-symmetric about the ring midline. Then 



lN+j 



IN 



-{I 



N- 



In), 



(12) 



and if j < N and k > N, then 
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j — N N — k 



(13) 



It is straightforward to show that by virtue of identities (11)— (13), the summation in equation 
(9) can be rewritten in terms of the {xi} as 



2N-1 , 



i=j—l l=i 



' k-j 



= E x *E 



i=i /= 

i=N-l 



(14) 



+ E 



2N-1 , 2=V , 

7k~7V 2^ (1 _ 



i=AT-l 

+ E ^E 



f(2N-j-iy 



When equation (14) is substituted into equation (9), we obtain a linear system of iV — 1 
independent equations for N — 1 unknowns, {xj}. These are solved numerically using the 
LAPACK subroutines. The set of inclinations, {/,■}, are derived from {xj} via (10). 

Figures 5 and 6 summarize the results of our calculation for the a and (3 rings of Uranus. 
These are the only narrow rings whose mean inclinations, J, have been measured. We predict 
inclination profiles that are much flatter than those predicted by the standard self-gravity 
model. The larger ring masses that arise from the inclusion of interparticle collisions enable 
rings to maintain nodal alignment by warping their geometry by relatively small amounts. 
The peak-to-peak amplitude of the warp, defined as a(l2N-i—h), is ~7 meters for the a ring 
and ~3 meters for the j3 ring. By contrast, the standard self-gravity model, which employs 
wire masses derived without regard to interparticle collisions, predicts warp amplitudes of 
order 1 km. 

Figures 7 and 8 describe inclination profiles for the Maxwell and Colombo ringlets of 
Saturn. The mean inclinations, J, of these rings with respect to the local Laplacian plane have 
not been measured. Fortunately, we have found by numerical experiment that the fractional 
variation of inclination across the ring, (I — I)/ 1, is independent of the mean inclination. 
The top panels of Figures 7 and 8 should therefore be fairly accurate. Numerical experiment 
also reveals that the absolute value of the warp across the ring, a(l2N-i — h), scales linearly 
with the unknown /; in computing the absolute value of the warp for the bottom panels, we 
have simply guessed I = 3 x 10~ 4 rad for both rings, values comparable to their observed 
mean eccentricities. 
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Fig. 5. — Inclination profiles of the Uranian a ring using two different mass models. Right- 
hand panels portray the inclination profile using the surface density profile computed in 
§2.1, which accounts for ring self-gravity, the planetary quadrupole field, and interparticle 
collisions. Left-hand panels correspond to inclination profiles computed under the standard 
self-gravity model that does not account for interparticle collisions. Upper panels describe 
the fractional variation in inclination across the ring, while lower panels plot the amplitude 
of the vertical warp across the ring. The large masses predicted by our work yield much 
flatter rings than those predicted by the standard self-gravity model; we predict vertical 
warps of ~7 meters across a ring that is several kilometers wide. 
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Fig. 6. — Same as Figure 5, but for the (3 ring of Uranus. Warp amplitudes of only ~3 
meters are predicted by our models. 
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Fig. 7. — Same as Figure 5, but for the Maxwell ring of Saturn. While the mean inclination, 
I, is unknown for the Maxwell ring, its value does not affect our calculation of the fractional 
variation of inclination across the ring, (I — I)/ 1, shown in the upper panels. By contrast, 
the absolute value of the vertical warp, a8I = a(I — I), scales linearly with 7; to compute 
the values displayed in the lower panels, we have assumed I = 3 x 1CT 4 rad. Accounting for 
interparticle collisions raises the mass of the ring and therefore reduces vertical warping; the 
effect is less dramatic for the Maxwell ring than for the Uranian a ring because the radial 
width of the former is greater than that of the latter, so that the solution is less sensitive to 
boundary conditions. 
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Fig. 8. — Same as Figure 7, but for the Colombo (Titan) ring of Saturn. In generating the 
bottom panels for aSI = a(I — I), we have assumed I — 3 x 1CT 4 rad. The top panels 
displaying the fractional variation of inclination are independent of this assumed value. 
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Note that vertical warping is substantially greater for the Saturnian Maxwell and 
Colombo rings than for the Uranian a and (3 rings. Accounting for interparticle collisions 
near ring edges always increases ring masses and flattens ring warps, but such effects become 
less important as the ring width increases. The Maxwell and Titan rings are each more than 
25 km across at quadrature, while their Uranian counterparts are less than 7 km across. 
Collisional stresses near ring edges are not as effectively communicated across wider rings. 



We check a posteriori the validity of our approximation in §2.2 that interparticle colli- 
sions play no direct role in determining the inclination profile. Collisions matter most near 
ring edges (CGOO); an upper limit for the collisional acceleration felt by a ring particle lying 
with a mean free path, A ~ Q>/n, of the ring edge is c\j\ ~ Q>n. The component of this 
acceleration that is perpendicular to the orbit plane of this last ring particle is c^nqi, where 
qi = a(dl/da) is the inclination gradient evaluated at the edge. This normal acceleration 
generates a rate of nodal precession (relative to the collision-induced rate on the ring midline, 
which is zero) of ^CbTiqi / rial = Cbqi/al. 

We divide this maximum, collision-induced differential precession rate by the differential 
rate of this last endwire due to the quadrupole field of the planet. For collisions to be 
unimportant, this dimensionless number should be less than unity. For the Uranian a ring, 
it is 



where the value for qi at the ring edge is obtained from our computed inclination profile. For 
the (3 ring, we derive Co 0.2(c b /2 cms -1 ); for both the Saturnian Maxwell and Colombo 
(Titan) ringlets, Co ps 0.1(c b /2 cms -1 ). Note that these values of Co for the Saturnian 
ringlets are independent of the assumed / since qj/I oc SI/ 1 is independent of I. We conclude 
based on these smallish values of Co that collisions are of marginal direct importance to the 
inclination profiles of our ring models. As we recede from the ring edge towards the ring 
interior, collisions become even less important as surface density gradients decline. 



2.3. Neglect of Collisions in Vertical Structure 




(15) 
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3. STABILITY 

Here we prove that circular, nodally locked ringlets are linearly stable to perturbations 
to their inclinations and nodes. A corresponding proof for the stability of apsidally locked, 
co-planar rings would read similarly. We account not only for ring self-gravity and the 
planetary quadrupole field, but also for interparticle collisions, albeit in a simplistic way. 

We restrict ourselves to circular rings for simplicity and for this proof only. Our hope 
is that accounting simultaneously for ring eccentricity would not change our conclusions. 
While this hope must be tempered by considerations of a time-varying uo that we have not 
explored in depth, we take some consolation from the Appendix in which we show how a 
variable uj impacts force balance in the vertical direction by seemingly modest amounts. 

In §3.1, we derive the secular changes to the inclination and node of a test particle due 
to the gravitational influence of a nearby massive wire. In §3.2, we show how narrow rings 
respond to small perturbations by oscillating in a number of normal harmonic modes. 



3.1. One Wire and a Particle 

Consider the perturbations induced by a massive, circular wire on a nearby test particle. 
The equator plane of the planet defines the reference plane, and z is the unit vector perpen- 
dicular to this plane. Let a p , I p , and Q p be the particle's semi-major axis, inclination, and 
longitude of ascending node, while a, I, and Q refer to the corresponding orbital elements 
of the perturbing wire. Let $ denote the particle's angular position away from the node, 
$ = — Q p , where 9 is the particle's longitude. Define Aa = a — a p , AI = I — I p , and 
Af2 = f2 — Figure 9 illustrates the geometry. 

Assume that |Aa| <C a, |AJ| <C /, and |AQ| <C 1. If the wire has total mass m, then 
its linear density is 

p = m/27ia . (16) 

Because the separation between the particle and the wire is small, the wire can be approxi- 
mated locally as straight. The gravitational acceleration of the particle induced by the wire 
is therefore 

F = 2Gpd/d 2 , (17) 

where d is the perpendicular vector from the particle to the wire, d = |d|, and G is the 
gravitational constant. 

As depicted in Figure 10, the vertical distance from the particle to the wire at fixed 



-20- 




Fig. 9. — The three-dimensional geometry of the problem. The central planet rests at the 
origin, and the z-axis is aligned with the planet's spin axis. A test particle is located at 
point p and is displaced from its node on the planet's equator plane by an angle d is 
the perpendicular vector from the particle to the wire. Angles Q and / are the longitude 
of ascending node and inclination of the wire, respectively, Q p and I p are the corresponding 
angles for the particle, and AQ = Q — Q p . The inclinations, particle-wire separation, and 
misalignment of nodes have all been exaggerated for clarity. 
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Fig. 10. — A two-dimensional projection of Figure 9 that plots height above the planet's equa- 
tor plane versus longitude on that plane. To achieve an exaggerated difference in orbits, we 
make Af2 large and negative for this figure. Because the angle a is 0(1), d z = Az [1 — 0(I 2 )\. 
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longitude is 

Az = asm Jsin^ — Afl) — a p sin 7 P sin $ 

= A(ai")sintf-a/AQcostf, (18) 

where A(al) = al — a p I p . Borderies et al. (1983b) assume A(al) aA7, which is equivalent 
to Aa/ a <C A// J; since we have derived inclination profiles in §2.2 for which their assumption 
is not valid, we keep our more accurate expression. Referring again to Figure 10, we note 
that since 0(a) = 0(1), 



d z = d • z = Az[l - 0(I 2 )} . (19) 

Furthermore, 



d 2 = d\ + Aa 2 . (20) 

The changes in I p and Q p due to F are given to leading order by (see, e.g., Murray & 
Dermott 1999) 

^E(tf) = ^ (21a) 

^(tf) = ^2!* (21b) 
dt n p a p 

where n p is the particle's mean motion, F z = F • z = 2Gpd z /d 2 ~ 2GpAz/d 2 is the vertical 

component of the perturbing force, and we have neglected terms of order I 2 . Substituting 
eqs. (16), (17), (18), and (20) into (21a) and (21b) yields 

dn v , x ^sin 2 tf- Afisintfcostf , x 

— ^(0) = C — ^ 5 , (22a) 

d * i+f^lsin^-^cos^ 

V Aa Aa J 

dl v , ^sin^cos^-JAficos 2 ^ , s 

-TrW = C a — -3, (22b) 



where 



l+(^sin^-^costf) 



m / a \ 2 



C= — — — ra (23) 

and n is the mean motion at semi-major axis a. Where appropriate, we have neglected 
differences between a p and a. 
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To obtain the secular changes in I p and Q p , we average dQ p /dt and dl p /dt over one full 

orbit of the test particle; to wit, X = (l/2n) X(-&) d-d. To simplify this integral, first 
define Q and 7 such that 

Qcos 7 = (24a) 

^ . al AQ 

Q Sm 7 = ~Aa~' ( } 
Then eqs. (22a) and (22b) can be re-written 

dn P(A\ _ n( Aa \ Qsin^sin(^-7) 
_ nf Aa \ Qcos^sin(^-7) 

" l + Q 2 sin 2 (tf- 7 )' ( } 

Averaging these equations over 1? yields 

^ = C^V(-g 2 )gcos 7 = ^(-Q 2 )^ + ^ (26a) 



dt J G \ al J \ I 

*Iz) = -C (— W-Q 2 )Q sin 7 = -CH(-Q 2 )IAQ, (26b) 



dt J G \ a 



where the subscript G denotes wire gravity and 



)2 , _ v'<F~T 1 



H(-Q 2 ) = v " ' (27) 
The forms of C and if are chosen to be consistent with eq. (6) of Borderies et al. (1983b). 



3.2. N Self-Gravitating Wires with Oblate Planet and Collisions 

Now consider the interaction between N circular wires. We demonstrate that the equi- 
librium configuration of the wires is linearly stable. For this purpose, we do not take the 
wires to be equally spaced in semi-major axis, but rather define their spacing such that all 
wires have equal mass; the wires must therefore be more closely packed in regions of higher 
surface density. This de-composition by mass rather than by semi-major axis is critical to our 
proof of stability. Denote the total ring mass by M, so that each wire has mass M/N. Let 
wire j have semi- major axis a,j, longitude of ascending node Qj, and inclination Ij. Define 
mean elements by a = aj/N, Cl = £lj/N, and I = J] Ij/N; let n denote the Keplerian 
mean motion at semi-major axis a. As before, we assume that variations in orbital elements 
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across the ring are small: for any wire j, \dj — a\ <^a, \Ij — I\ <^ I, and \Qj -Q| C 1. We 
adapt eqs. (26a) and (26b) to write down the time rates of change of £lj and Ij due to the 
gravitational attraction of wire k: 

dtlA = C Jk H(-Q%)(^ + ^) (28a) 



^ ' G(k) 



dt 



G{k) 



-C jk H(-Q 2 jk )I(Q k - Qj) , (28b) 



where 



C lk = C kj = M ( — - — ^ n (29) 
3 3 TiMpN \a k — dj J K ' 

a k I k - ajlj \ 2 / aI(Q k - Qj) x 2 



Q% — Qkj — ( J J + V K " ) ■ (30) 

J J \ a k - aj J \ a k - cij J 

Interparticle collisions are assimilated as follows. According to equation (2), the col- 
lisional repulsion between two neighboring wires is inversely proportional to the distance 
between them, and the resulting acceleration is directed oppositely to that of inter-wire 
gravity. This prescription implies that collisional forces act effectively as anti-gravitational 
forces, having the same dependence on inter-wire separation as self-gravity but with an op- 
posite sign. Thus, we can account for collisions by prepending a factor gj k to the right-hand 
sides of eqs. (28a) and (28b). The factor g jk represents the effects of collisions; the particular 
value of Qj k depends on the local velocity dispersion; it is always < 1, and it is negative 
where collisional repulsion is stronger than gravitational attraction. 

Finally, consider the central planet's quadrupole field, which causes no secular change 
in {Ij} but induces precession of the nodes {£lj}. The effect is described to leading order by 

(1r) = -\^i^M-^ 712 ■ (3D 



A linear expansion yields 



— 2 

dQj\ dQ 21 i R n \ _ o>i — a 



T + T^P ^ (32) 



dt J q dt I 

where dCl/dt = -\J 2 R 2 v ^fGM P ~ar 1 1 2 . 

The total rates of change of flj and Ij due to all other wires read as 
dQi 



dt 
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dVt 21 f R«\ 2 a- — a , 
+-^ + -rJA^)n^— (33a) 



dt 4 

^ = -^^C, fe //(-gy/(fi fe -Q,). (33b) 



Eq. (33b) implies that / is conserved. Now denote the equilibrium inclination of wire j by 
Ij, and define 

/ - 1° 

*i = v • (34) 

Equations (33a) and (33b) are expressed more concisely in terms of 

^ = ^g jk C jk H(-Cf jk )(* k -* j ) + ^ (35a) 
^ = -J29 jk C jk H(-Q%)(n k -^). (35b) 

k^j 



The terms in eq. (33a) proportional to (a k — dj)/a and (a^ — a) /a have been absorbed into 
Ij (i-e., $,) and 4° (i.e., $ fe ). 

As a final simplification, define iV x 1 column vectors x, y, and u with elements = Qj, 
Uj = $j, and ttj = 1. Equations (35a) and (35b) can then be written in matrix form: 

d - v~ , \ 

— x = iiyH — —u, (3oa) 

|y = -Xf, (36b) 
where K is an N x N symmetric matrix whose off-diagonal entries are given by 

Kjk = K kj = g jk C jk H(-Q 2 jk ), (37) 
and whose diagonal entries equal 

K 30 = -Y^mkC 3k H{-Q) k ). (38) 

Note that K is symmetric only because we have chosen our wires to have equal mass. 
Furthermore, we can take H to be constant because the fractional variation in this quantity 
is of order Q 2 times the fractional variation in Q, and Q ~ C 1 from our equilbrium 
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model. Finally, we assume that variations in g^ are small. Then K is a constant matrix 
and eqs. (36a) and (36b) are linear. 

Since K is symmetric, it possesses N orthogonal eigenvectors v , . . . , vn-i with real 
associated eigenvalues Xo? • • • , Xn-i- Inspection reveals that v = u and Xo — 0. Thus the 
complete solution to eqs. (36a) and (36b) is 

N-l ^ 

n=l 
N-l 

V = ^ c n {7„ cos [x„(* - 0n)] , (39b) 
Jl=l 

where the 2A^ — 2 constants Ci, . . . , ca?_i and 0i, . . . , 0jv-i are determined from the initial 
relative inclinations and longitudes of ascending nodes. The two remaining constants are to 
and 7, where to is set by the mean longitude of ascending node. Equations (39a) and (39b) 
imply that an inclined, circular ring responds to small perturbations to the inclinations and 
nodes of its constituent ringlets in the presence of self-gravity, planetary oblateness, and (a 
simplified prescription of) interparticle collisions by oscillating in a number of normal modes. 

These modes resemble the vibrations of a string with free boundary conditions and 
variable density. A preliminary numerical study in which we modelled the ring with = 1 
suggests that the fundamental (lowest, non-zero) frequency is 

where Aa is the full width of the ring. For this mode, half the ring precesses nodally in one 
direction relative to the ring midline, while the other half precesses in the other direction. 
For the Maxwell ringlet, 2ft /xi ~ 8 yr. This same study indicates that Xj — ^Xij for j > 10. 
We will use this numerical model in our discussion below regarding the possibility of exciting 
a high order mode by using force variations that arise from a time-varying u (see also the 
Appendix). 



4. SUMMARY AND DISCUSSION 



Our crude and prescriptive treatment of streamline dynamics near resonantly perturbed 
ring boundaries prevents us from claiming great accuracy in our computed surface density 
and inclination profiles. Nonetheless, we believe the following conclusions to be robust: 
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1. Narrow rings have total masses of order 10 19 g. 

2. Their surface densities increase towards their edges. 

3. The Uranian a and f3 rings are vertically warped by a height of order 10 m. The 
radially wider Saturnian Maxwell and Titan rings are warped by a height of order 1 
(7/3 x 10~ 4 rad) km. An equivalent and more model-independent statement is that the 
fractional variation of inclination is order 10~ 3 across the Uranian ringlets and of order 
1CT 1 across their Saturnian analogues. 

4. Nodally (apsidally) locked rings are linearly stable to perturbations to their inclinations 
(eccentricities) and nodes (apses). 

Analysis of stellar occultation data by the Cassini spacecraft may test predictions 1-3 for 
the Maxwell and Colombo ringlets of Saturn. 

For a ring to maintain apsidal and nodal alignment, its argument of pericenter must 
precess. This means that where the ring is "pinched" — in other words, its pericenter, the 
narrowest portion of the ring, where self-gravitational and collisional forces are strongest— 
circulates relative to the node of the ring on the equator plane of the planet. Such circulation 
must introduce time-variable forces into the vertical equations of motion for the ring; these 
time-variable forces have a frequency of 2uj, or about 0.08 cycles per day for the Saturn 
Maxwell ringlet. Might these time- variable forces drive bending waves across the ring? First 
note that the driving frequency is more than 2 orders of magnitude greater the fundamental 
bending frequency, Xi- Therefore any mode that is excited will have a wavelength consider- 
ably shorter than the overall width of the ring. A preliminary study suggests that the radial 
wavelength of the bending mode whose natural frequency matches the driving frequency is 
of order one kilometer for the Saturnian ringlets, a length scale possibly within reach of 
the Cassini camera. We defer a more thorough investigation of this possibility, including 
estimates of the amplitude of the mode, to future study. In the event that Cassini finds the 
vertical structure of a narrow ring to vary over a period of a few weeks and over kilometer 
lengthscales, one might look to the periodic stresses associated with a time-varying uj for 
the cause. Ultimately, measurements and a theoretical understanding of any modes excited 
might constrain the mass of the ring independently of optical depth soundings. 

We have solved for the inclination gradient of a given ring by first finding the surface 
density profile. A key input for this procedure was the eccentricity profile of the ring, 
e(a), which we obtained from the (spatially averaged) observations. While this procedure 
is physically self-consistent, it begs the question of where the eccentricity profile originated. 
Another way of saying this is to note that our equilibrium solutions are not the only ones 
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possible in principle. For example, the steps of our procedure could be reversed: we could 
begin by positing an inclination gradient, solve next for the masses required to maintain 
nodal alignment, and conclude by deriving the eccentricity gradient required to maintain 
apsidal alignment. We did not adopt this procedure because the observations supplied 
the eccentricity gradient and not the inclination gradient. Why narrow rings self-organize 
themselves to exhibit dimensionless eccentricity gradients of order q e = ade/da ~ 0.5 and 
dimensionless inclination gradients that are significantly smaller — of order qj = adl /da ~ 
lO^-lCT 3 for the Saturnian cases and of order 10" 3 -10~ 5 for the Uranian cases — is unclear. 

Many issues remain unresolved in our understanding of narrow rings. In the more than 
twenty years since their discovery, we still ask the following questions: 

1. Are narrow rings primordial? According to classical theories of ring shepherding, at- 
tendant satellites extend the viscous spreading time of a ring by a factor of order the 
mass ratio between the shepherd and the ring. This extension factor is modest for 
the Uranian e ring, whose mass is thought to be comparable to that of its shepherds 
(CG00). The situation is likely to be similar for other rings, especially given the large 
ring masses that we are deriving. Estimates for the spreading time of the e ring range 
from 3 x 10 4 yr to 2 x 10 8 yr, all of which are shorter than the age of the solar system 
(Chiang 2003). But if narrow rings are not formed from primordial circumplanetary 
disks, how can we explain the tiny mean inclinations exhibited by the Uranian rings 
with respect to the equator plane of their severely oblique host planet? New estimates 
of ring lifetimes based on non-classical theories of shepherding (see, e.g., Goldreich & 
Porco 1987, and Borderies, Goldreich, & Tremaine 1984) would be welcome. 

2. Why do all narrow rings exhibit positive and not negative eccentricity gradients? 

3. Why is the inclination gradient smaller in magnitude than the eccentricity gradient? 

We might hope to answer these questions by modelling more carefully the distortion of 
streamlines near ring edges by shepherd satellites. 

It is a pleasure to thank Joanne Cohn for introducing C.J.C. to E.I.C. We thank also 
Carolyn Porco for faxing pages from her thesis to us, and an anonymous referee for a helpful 
report that prompted us to examine our statements regarding bending modes more carefully. 
This work was supported in part by a Space Grant Summer Fellowship awarded to C.J.C. 
by the Space Sciences Laboratory at Berkeley. E.I.C. acknowledges support by National 
Science Foundation Planetary Astronomy Grant AST-0205892 and Hubble Space Telescope 
Theory Grant HST-AR-09514.01-A. 
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A. Time-varying uj 

Here we explore the effect of a time- varying argument of pericenter, u, on our equation 
of motion (6) for nodal precession in the presence of ring self-gravity. Consider two elliptical 
wires orbiting a spherical planet, one placed in an arbitrary reference plane, and a second 
placed at an inclination, Al, relative to that plane. Take the longitudes of pericenter of the 
two wires to be equal, uj\ = u; 2 = 0. The longitude of ascending node of the second wire is 
f2. How sensitive is the orbit-averaged rate of nodal precession of the second wire, (ti), to 

n = -uoi 

At every true anomaly, /, of the second wire, the perpendicular distance between the 
second wire and the first is 



y(Aa) 2 (l -g e cos/) 2 + a 2 (A/) 2 sin 2 (/-fi), (Al) 

where we have assumed that the difference in semi-major axes of the wires, Aa, is much 
smaller than the semi-major axis of either wire, a. The eccentricity gradient between the 
wires is q e = aAe/Aa, where Ae is the difference in wire eccentricities. The first term under 
the square root in equation (Al) represents the in-plane separation between the wires, while 
the second term represents the out-of-plane contribution. Re- write (Al) as 

d » Aa^j (1 - q e cos f f + qj sin 2 (/ - Q) , (A2) 

where qj = aAI/Aa. 

By Gauss's equation, the instantaneous rate of nodal precession of a test particle whose 
orbit coincides with the second wire reads 

U " naAI ' {A6) 

where n equals the mean motion of the second wire, F z equals the vertical acceleration felt 
by the test particle, 

2GpaA/sin(/-^) 
F z = j2 , (A4) 



and p equals the linear mass density of the first wire of mass m, 
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We have dropped corrections to (A5) due to wire eccentricity because the acceleration (A4) 
is already small in AI. Time-averaging Cl over one orbital period yields 

to = -^hh B{ *^> (A6) 

b=- r ^v-v) df (A7) 

2nJ (l-g e cos/) 2 + g 2 sin 2 (/-fi) 1 ' 

Equation (A6) is identical in form to equation (6). The sensitivity of (Cl) to w = — fl is 
contained in the integral, B. Numerical integration of (A7) reveals that B e [0.613,0.896] 
for q e = 0.5 and qj = 0.1. If q e = 0.5 and qi = 0, then B G [0.619,0.920]. The variation 
of force with varying uo is restricted to less than a factor of 1.5. In §2, we fixed B = 0.77 
for a maximum 24% fractional error; the error so accrued is less than the error introduced 
by our prescriptions for force balance near ring boundaries. However, the force variation is 
large enough that it may excite observable short-wavelength bending waves across the ring, 
as we discuss in §4. 
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Table 1. Observed and Derived Parameters of Narrow Eccentric Rings 



Ring Planet a(km) a Aa(km) a e(xl0 3 ) a (g e ) a 7(xl0 3 ) a (g/) b (AI/I) b M(10 19 g) b 



a 


Uranus 


44718 


7.15 


0.761 


0.472 


0.265 


0.00095 


0.00057 


0.467 





Uranus 


45661 


8.15 


0.442 


0.370 


0.089 


0.00037 


0.00075 


0.495 


Maxwell 


Saturn 


87491 


64 


0.34 


0.46 


(0.3) 


(0.070) 


0.17 


5.00 


Colombo (Titan) 


Saturn 


77871 


25 


0.26 


0.44 


(0.3) 


(0.055) 


0.059 


1.64 



a Uranian ring values taken from Tables I and VII of French et al. (1991). Saturnian ring values taken from Table Hid 
of Porco (1983) and Table 1 of Porco (1990). Values for T in parentheses are guessed. 

b Derived from this work. Values in parentheses for qj = aAI/Aa, where AI is the difference in orbital inclinations 
between inner and outer ring edges, scale linearly with the guessed value of /. 



Table 2. Theoretical Model Parameters 



Ring 


Cj (cms x ) 


Q>(cms r ) 


w r (km) a 


A(km) b 


N 


a 


0.1 


2.0 


0.44 


0.079 


1000 


P 


0.1 


2.0 


0.45 


0.081 


1000 


Maxwell 


0.1 


3.0 


0.56 


0.126 


3000 


Colombo (Titan) 


0.1 


2.0 


0.50 


0.071 


3000 



a Resonant width, computed as w r = ayj M sat /M P , where the mass of an 
individual (as yet unobserved) shepherd satellite equals M sat = 8.4 x 10 18 g 
for the Uranian ringlets and M sat = 2.3 x 10 19 g for the Saturnian ringlets. 

b Mean free path near ring edge, computed as A = Q,/n, where n is the 
average mean motion of the ring. 



